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Abstract. - We study the relaxation modes of an interface between a lyotropic lamellar phase 
and a gas or a simple liquid. The response is found to be qualitatively different from those of 
both simple liquids and single-component smectic-A liquid crystals. At low rates it is governed 
by a non-inertial, diffusive mode whose decay rate increases quadratically with wavenumber, 
— Aq 2 . The coefficient A depends on the restoring forces of surface tension, compressibility 
and bending, while the dissipation is dominated by the so-called slip mechanism, i.e., relative 
motion of the two components of the phase parallel to the lamellae. This surface mode has a 
large penetration depth which, for sterically stabilised phases, is of order (dq 2 )' 1 , where d is 
the microscopic lamellar spacing. 



Lyotropic lamellar phases occupy large portions in the phase diagrams of amphiphilic 
molecules (surfactants) in solution [1,2]. They consist of stacks of parallel fluid membranes 
separated by microscopic (~1-10 nm thick) layers of solvent [fig. ^a)], thus having the sym- 
metry of a smectic-A liquid crystal [3]. These phases appear in numerous applications, e.g., 
in the cosmetic and detergent industries. Lamellar bodies are found also in biological systems 
such as the lung [4]. Apart from the free surfaces of lamellar phases with air, surfactant 
phase diagrams contain also large coexistence regions (so-called immiscibility gaps) [1], in 
which a lamellar phase has an equilibrium interface with an isotropic liquid such as a dilute 
micellar solution or an L3 (sponge) phase [1,2]. Spherical lamellar structures in the form of 
multilamellar vesicles (onions) dispersed in a solvent are commonly encountered as well [5]. 
Thus, the surface response of these phases is a fundamental issue relevant to a large variety of 
experimental systems. The static response of smectics to surface deformations was thoroughly 
studied [3,6,7]. The dynamics of surface perturbations in thermotropic (single-component) 
smectics were investigated as well, for both semi-infinite systems and finite films [8-16]. In 
the current Letter we analyse the surface relaxation of lyotropic (two-component) lamellar 
phases and demonstrate the essentially different surface dynamics of this ubiquitous class of 
materials. 

Surface modes characterise the relaxation dynamics of surface perturbations whose am- 
plitude decays with increasing distance into the bulk material [fig-^a)]. In a simple liquid, 
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Fig. 1 - a) Schematic view of the system. A semi-infinite lamellar phase occupies the region z < 0, 
having an interface at z — with a gas or a simple liquid. A surface perturbation is shown, whose 
amplitude decays to zero at z — > — oo. b) Visualisation of the slip flow, in which the two components 
of the lamellar phase have different average velocities. 



having mass density p, viscosity rj and surface tension 7, the surface dynamics depend on two 
dimensionlcss numbers: the Reynolds number, Re = p\ui\ / (rjq 2 ) , and the Capillary number, 
Ca = 7/|a;|/(7|g|), where q and u) are the perturbation wavenumber and frequency, respec- 
tively. If inertia is neglected, one is left with the Capillary number alone, leading to a decay 
rate ~ j\q\/r]. However, substituting this result back in the requirement Re <C 1 yields 
\q\ 3> pj/r/ 2 . Thus, e.g., for water such a non-inertial, overdamped response is restricted to 
wavelengths much smaller than 1CP 2 pm, while larger-wavelength perturbations are governed 
by inertial capillary waves [17], whose dispersion relation is set by ReCa ~ 1. 

The viscoelastic hydrodynamics of smectic-A liquid crystals [3] were formulated in the 
seminal work of ref. [18]. Due to the anisotropy of these materials the viscosity tensor contains 
five independent coefficients which, for incompressible flows, can be reduced to three, denoted 
by ?7m, Vt and 7/v [19]. While tjm is related to relative sliding of layers, and is therefore 
comparable to the viscosity of the solvent (water) in lyotropic phases, the other two are 
associated with distortions of the membranes themselves and, hence, are comparable to the 
overall apparent viscosity of the liquid, which is 2-3 orders of magnitude larger [19, 20]. The 
elasticity is characterised by a compression modulus B and a bending modulus K. These 
elastic moduli define two Deborah numbers in addition to the Capillary and Reynolds numbers. 
The dynamics thus depend on four dimcnsionless numbers, which we define as 

De B = 7]m\u)\/B, Dc k = Ji M \uj\/(Kq 2 ), Ca = ^mM/^M), Rc = p\ui\ / (r] M q 2 ) . (1) 

(More accurately, in the case of lyotropic phases one should consider separately, depend- 
ing on the frequency, two layer-compression moduli — the modulus at constant surfactant 
concentration (for high frequencies) and the modulus at constant surfactant chemical po- 
tential (low frequencies), which is typically much smaller [21]. Our main focus here is on 
the slow response for which B should be regarded, unless otherwise noted, as the smaller, 
fixed-chemical-potential modulus.) 

The surface modes of a semi- infinite thermotropic smectic phase were studied in ref. [14]. 
Despite the much higher viscosity of these phases compared to ordinary simple liquids (water) , 
the conclusion of this work was that the low-frequency surface modes were inertial, elastic 
(Rayleigh) waves of second sound, their dispersion relation being set by ReDcs ~ 1. In a 
detailed analysis of low-Reynolds-number surface modes of thermotropic smectics [22], we 
have reached a similar conclusion, namely, that non-inertial relaxation is restricted to small 
(sub-micron) wavelengths. On the other hand, it is hard to envisage the surface dynamics 
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of a small lamellar droplet or an onion (of, say, 10-10 2 /im scale), whose inter-layer spacings 
and surroundings are pervaded by viscous solvent, as dominated by underdamped inertial 
waves. The disagreement between this intuition and the aforementioned results suggests that 
lyotropic lamellar phases may be dominated by another dissipative mechanism, different from 
the usual viscous one. 

The hydrodynamic theory of smectics was extended to lyotropic phases in refs. [19,23,24]. 
Brochard and de Gennes [19] were the first to recognise that, unlike isotropic binary mixtures, 
the spatial organisation of the two components in a lamellar phase allows for collective motion 
of one component relative to the other parallel to the layers [fig-IHb)]. This bulk slip mode 
was observed experimentally [23,25-28]. We show in the current work that the additional 
dissipation due to slip has profound consequences for the surface dynamics. The dominance 
of this dissipative mechanism can be demonstrated by the following simple argument (to be 
made rigourous later on). The transport coefficient associated with slip, /i, relates to the 
mobility of a Poiseuille flow of solvent between two membranes [fig-^b)], /i ~ d 2 /rjw> where 
d is the inter-membrane spacing and rj^ the solvent viscosity [19]. The resulting friction force 
per unit volume is ~ H~ 1 v, v being a characteristic relative flow velocity, while the force 
density arising from viscous stresses is ~ ijryq 2 v. Hence, the slip dissipation dominates for 
q <C (t7t,vm) -1 ^ 2 ~ (lOd)" 1 — a condition satisfied for all physically relevant wavelengths. 

The system under consideration is schematically depicted in fig. 1(a). A semi- infinite 
lamellar phase occupies the region z < 0, the layers lying, on average, parallel to the xy plane. 
At z = there is an interface between the lamellar phase and a dilute gas or a simple liquid. 
We use Latin indices to denote vector components along all three axes (x, y, z) and Greek 
indices (a, /3) for components along the two lateral directions (x, y). 

The hydrodynamic description of thermotropic smectics requires six scalar hydrodynamic 
fields [3, 18]. Assuming an isothermal, incompressible flow, one is left with only three indepen- 
dent fields, e.g., two of the three components of the flow velocity i>j(r, t) at position r and time 
t, and an additional scalar field, u(r, t), defining the local displacement of the layers from their 
equilibrium position. Three dynamic modes result: two second-sound waves and a transverse 
shear (vorticity) diffusive mode. In lyotropics another scalar field is required [19,23], e.g., 
the local surfactant volume fraction. As a result, another mode emerges, i.e., the slip (baro- 
clinic) mode. In writing the hydrodynamic equations we use the variables Vi(r,t), u(r,t), and 
a pressure field p(r,t) (of which, we recall, only three are independent). To include the slip 
degree of freedom in lyotropics another velocity field is added [19], which we take as the lateral 
velocity of the surfactant component, v s _ a (r,t). Velocity differences between the two compo- 
nents, v s , a — v a , are driven by in-plane stresses in the membranes [19]. We focus here on the 
simplest treatment, where the membranes are considered as incompressible two-dimensional 
(2D) liquids with no 2D viscous stresses. (This assumption is also consistent with the limit, 
mentioned above, of the fixed-chemical-potential compression modulus being much smaller 
than the fixed-concentration one.) In this description the driving force for slip is provided 
solely by a lateral pressure gradient, d a ip(r,t). Out of the three new variables, i> SiQ and if), 
similarly, only one is independent. 

In terms of the aforementioned variables, the hydrodynamic equations for a lamellar phase 
take the form 



p{d t Vi + vjdjVi) = djCTij 
diVi = 
d a v s . a = 

V s , a -V a = -fid a 1p 



(2) 
(3) 
(4) 
(5) 
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d t u -v z = 0. (6) 

Equation J2J i s the Navier-Stokes equation, where try, to be specified below, is the stress 
tensor. Equations Q and J1J impose mass conservation on the incompressible overall flow and 
2D surfactant flow, respectively. Equation JSJ asserts a linear relation between the slip velocity 
and its driving force. In eq. 10 we have assumed that permeation of material perpendicular 
to the layers is negligible [3]. We write the stress tensor as a = a 1 ' + a d , where the "reactive" 
part is given by [19,29] 

o"L = a \y =-p-tp + Kd,3p z u, o\ z = -p + Bd z u, a r az = -Kdpp a u, a T xy = 0, (7) 

and the dissipative part by [19] 

a £p = VridaVp + dpv a ), <J ZZ = 2i]yd z v z , a A az = r/ M (d a v z + d z v a ). (8) 

At z — the tangent force per unit area, a az , must change continuously across the 
interface, whereas the perpendicular one, a ZZl has a discontinuity due to surface tension. If 
the surface is free (the other phase is a dilute gas), these conditions reduce to 

z = : <7 az = 0, a zz = ~fd aa u. (9) 

Since we deal with surface perturbations, the bulk far away from the interface is assumed to 
remain at equilibrium, 

Z — > — OO : Vi — v s ,a = u = 0. (10) 

If the phase in the region z > is a liquid, its hydrodynamic equations should be added, and 
the boundary conditions at z = are modified. Since the introduction of a simple liquid leads 
to a minor effect [22] , we shall briefly comment on it below and otherwise focus on the simpler 
case of a free surface. 

We are concerned with small deviations from equilibrium and, therefore, the nonlinear term 
in eq. © can be neglected. Equations l(2jl-l(TU)l then define a closed set of linear equations 
with the appropriate boundary conditions, from which the surface modes can be derived. For 
simplicity we restrict the discussion to variations in the x and z directions only. Substituting 
a perturbation of the form f(x, z, t) oc e lg!C_lwt , where / stands for any of the hydrodynamic 
variables, we obtain a 4th-order equation for v z , 

d zzzz v z -(S + Q-2- Be^ - Re)q 2 zz v z + (1 - De^ 1 - Re)q 4 v z = 0. (11) 

In eq. (|Tl"1) the dynamic numbers (Deg, Dejf , Re) are as defined in eq. Q with \co\ replaced 
by iu>, and S and are two additional dimensionless parameters, the former associated with 
slip, S = (77mW 2 ) -1 ! an d the latter with the viscosity anisotropy, 9 = 2(ryx + Vv)/vm ^> 1. 
The solution to eq. Ijllfl is given by 

v z = (C+e a+z +C-e a - z y qx - luJt , a± = \q\ [(S + Q-2 - Dc^ 1 - Re ± T 1/2 )/2] 1/2 
T = (5 + 9- 2- Dc^ 1 - Re) 2 - 4(1 - De^ 1 - Re), (12) 

where the two spatial decay coefficients, a±, have been chosen as the ones with a positive real 
part so as to satisfy the boundary conditions at z — > — oo, cq. (|10|l . Imposing the interfacial 
boundary conditions © yields two linear equations for the amplitudes C+ and C_, whose 
determinant is set to zero to find the dispersion relation, u> = u>(q). The resulting equation is 

Ca _1 (a + + a_)/|g| + Rc - (1 - Dc^ 1 - Re) 1/2 (S* + 9 - Dc^ 1 - Dc^ 1 - Re) = 0, (13) 
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where Ca is as defined in eq. Q with the replacement \lo\ — > iu). 

Equation (|13|l . together with the expressions for a± given in eq. (|12|l . can be solved for u> 
only numerically. However, we can significantly simplify the problem by noticing the following. 
First, the slip mechanism will dominate over the viscous one when S » 0. As already 
discussed above, this amounts to a condition for q which is practically always valid. Secondly, 
for typical lamellar phases the length A = (K/B) 1 ^ 2 is very small (of order a few d) and, as a 
result, De^ 1 ^> De^- 1 for any relevant q. We also take Re <C De^- 1 < 1, <C De^ 1 , and S > 
Dc^ 1 . (The self-consistency of these conditions will be later verified.) Under these assumptions 
we find a+/\q\ ~ (S - De^ 1 ) 1 / 2 > a_/|g| ~ [(1 - De^)/(5 - De^ 1 )] 1 / 2 . Consequently, eq. 
ifH?)) is simplified to 

Ca" 1 - [(1 - De^)(S - De^ 1 )] 1 / 2 = 0, (14) 

whose solution is 

oj = -iAq 2 , A = [Bfi + K/ m + y/(B[i - K/ m ) 2 + W/t7m]/2. (15) 

Equation l|15fl is our central result. It describes a decay rate of a diffusive mode, which 
is affected by all three restoring forces and the slip dissipative mechanism. To make these 
results more transparent, let us specialise to sterically stabilised lamellar phases, where B ~ 
(k B T) 2 /(nd 3 ), K ~ n/d, and 7 ~ k B T/d 2 , k B T being the thermal energy and k the bending 
modulus of a single membrane, typically a few k B T. (The low tension values in lyotropic 
phases, 7 ~ k B T/d 2 , arise from the diffusive contact of the surface membranes with membranes 
in the bulk, which are essentially tensionless.) Recalling that ryiviA* ~ d 2 , we find for such phases 
that all the terms appearing in eq. (fT5)l are of the same order and A ~ k B T / (riyid) . This yields 
10 2 -10 3 pm 2 /s for d — 10 nm and ?7m = 1CP 2 poise. The spatial decay coefficients in this case 
are a+ ~ d" 1 and a_ ~ dq 2 . Thus, the surface mode contains a boundary layer of microscopic 
thickness ~ d, and a deeply penetrating undulation whose thickness is much larger than the 
wavelength. The amplitude ratio of these two parts is C+/C_ ~ —a 2 _/q 2 : which is a small 
(negative) number, i.e., the dominant part is the deeply penetrating one. 

Upon substitution of the results in the initial assumptions, their validity is readily con- 
firmed. In particular, the Reynolds number is Re ~ pk B T / (rf^d) ~ 10 -4 -10 -3 for the above 
parameters and p = 1 g/cm 3 , irrespective of q. In the case of stiffer (electrostatically sta- 
bilised) phases where B ^> {k B T) 2 /(nd?), however, the assumption Re <C De^ 1 is satisfied 
only if B -C (K j p) 1 ! 2 j p. For typical parameters this implies B <C 10 9 erg/cm 3 . Hence, for 
very stiff lamellar phases inertia may become important, and one should return to the more 
general equations l(T^|l and (fH5)l . 

Static surface deformations of smectics have a large penetration depth as well, (Ag 2 ) -1 ~ 
(dq 2 )" 1 [6,7]. Thus, the static and dynamic penetrations turn out to be similarly deep despite 
their different physical origins. While in the static case the large depth arises from the fact 
that the layers are much more resistant to compression than to bending, in the dynamic case 
it stems from the high resistance to relative motion of the two components compared with the 
viscous resistance to their overall flow. Once a deeply penetrating undulation exists, a thin 
boundary layer must appear as well to satisfy the surface boundary condition for the transverse 
stress, <r az = 0. If the boundary layer were absent, the ratio between the two transverse-stress 
terms would be d a v z /d z v a ~ q 2 /a 2 ^ 3> 1, whereas in the presence of a boundary layer of 
thickness a^ 1 we have d a v z /d z v a ~ (g 2 /a 2 h )(C_/C + ) ~ q 4 /(a + a^) 2 ~ 1 ; thus allowing the 
transverse-stress contributions to balance one another. 

It should be noted that, having employed the assumptions of 2D-incompressible membranes 
and negligible Re, we have suppressed the second-sound and vorticity-diffusion modes, leaving 
only one bulk mode (the slip, baroclinic one) of the four existing in bulk lyotropic phases. 
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Rayleigh waves, therefore, cannot emerge from this analysis. Such modes should exist in 
lyotropics as in thermotropics [14]. Their frequency domain, however, is much higher than 
that of the relaxation discussed here. For the two rates to be comparable, one must have 
(B/py/ 2 q ~ Aq 2 , i.e., q^ 1 ~ (Bp) 1 / 2 /! ~ 1-10 nm. Hence, the time scales are well separated 
for any relevant q and the two responses can safely be studied as decoupled. (The separation 
of time scales is even larger, in fact, since one should use for the fast second-sound response 
the much larger compression modulus at constant concentration [21].) 

If the phase in the region z > is taken to be a simple liquid [22] , the only change is that 
the term CsT 1 in eqs. (|13[1 and (|14fl is replaced with (Ca" 1 — 2r]/rju), V being the viscosity of 
the simple liquid. Therefore, the simple liquid does not affect the above results so long as r\ <c 
7/(A|q|). For sterically stabilised lamellar phases, for example, this implies that only when 
the simple liquid is orders of magnitude more viscous than the solvent, r\ ~ ??m/(c%|) 3> tjm, 
will it have an appreciable effect on the surface relaxation. This is yet another consequence 
of the dominance of slip dissipation. 

To summarise, in contrast with simple liquids and thermotropic smectics, whose surface 
dynamics are governed by underdamped waves (capillary or Rayleigh waves, respectively), 
we have found that surfaces of lyotropic lamellar phases can relax via a much slower, over- 
damped diffusive mode over a wide range of wavelengths. The key ingredient underlying this 
qualitatively different behaviour is the presence of viscous solvent in between the membranes 
and the resulting relative slip of surfactant and solvent layers. The relaxation is remarkably 
slow due to the strong friction (small fi) introduced by this motion. For example, in sterically 
stabilised phases perturbations of a micron-scale wavelength are predicted to decay with a rate 
of ~ 10 3 s -1 . To our best knowledge, experiments concerning the surface modes of lamellar 
phases have not yet been performed. The predictions of this work should be verifiable in, e.g., 
dynamic scattering experiments. 

In this study we have considered semi-infinite, flat lamellar phases. Our results should be 
valid for finite films and curved lamellar structures {e.g., onions) provided that their size is 
sufficiently large, or the wavelength is sufficiently small. This restriction, however, is particu- 
larly severe in the current case because of the deep penetration of the inferred surface mode. 
The onion radius, for example, should be much larger than aZ 1 ~ (dq 2 )^ 1 3> q 1 - For a 
typical onion of radius R ~ 10 fim and d ~ 1 nm the wavelength must be much smaller than 
{Rd) 1 / 2 ~ 10 _1 fim. An extension of the theory to finite films and curved surfaces is needed, 
therefore, to accurately account for the surface dynamics of such lamellar objects. 

Finally, lamellar surfaces out of equilibrium have been long known to exhibit an intriguing 
instability involving multilayer finger-like structures (myelin figures) [30,31]. We hope that 
the elucidation of the relaxation of lamellar surfaces presented here will be instrumental also 
in the resolution of this long-standing puzzle. 
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